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Summary 

In the design and analysis of robust control systems for uncertain plants, representing the 
system transfer matrix in the form of what has come to be termed an “M-A model” has become 
widely accepted and applied in the robust control literature. The symbol M represents a transfer 
function matrix M(s) of the nominal closed-loop system, and A represents an uncertainty 
matrix acting on M(s). The nominal closed-loop system M(s) results from closing the feedback 
control system K (s) around a nominal plant interconnection structure P(s). The uncertainty 
can arise from various sources, such as structured uncertainty from parameter variations or 
multiple unstructured uncertainties from unmodeled dynamics and other neglected phenomena. 
In general, A is a block diagonal matrix, but for real parameter variations, A is a diagonal 
matrix of real elements. Conceptually, the M-A structure can always be formed for any linear 
interconnection of inputs, outputs, transfer functions, parameter variations, and perturbations. 
However, very little of the currently available literature addresses computational methods for 
obtaining this structure, and none of this literature (to the authors’ knowledge) addresses a 
general methodology for obtaining a minimal M-A model for a wide class of uncertainty, 
where the term minimal refers to the dimension of the A matrix. Since having a minimally 
dimensioned A matrix would improve the efficiency of structured singular value (or multivariable 
stability margin) computations, a method of obtaining a minimal M-A model would be useful. 
Hence, a method of obtaining the interconnection system P(s) is required. This paper presents 
(without proof) a generalized procedure for obtaining a minimal P- A structure for systems with 
real parameter variations. With this model, the minimal M-A model can then be easily obtained 
by closing the feedback loop. The procedure involves representing the system in a cascade-form 
state-space realization, determining the minimal uncertainty matrix A, and constructing the 
state-space representation of P(s). Three examples are presented to illustrate the procedure. 

1. Introduction 

Robust control theory for both analysis and design has been the subject of a vast amount 
of research. (See refs. 1 through 35.) In particular, robust stability and performance have been 
emphasized in much of this work, as, for example, in the development of H°° control theory. 
(See refs. 10 through 15 and 19 through 23.) Moreover, the development of robust control system 
design and analysis techniques for unstructured (refs. 1 through 9, 13, 19, and 21) as well as 
structured (refs. 16 through 35) plant uncertainty continues to be the subject of much research — 
particularly the latter. Unstructured plant uncertainty arises from unmodeled dynamics (such 
as actuator or engine dynamics) and other neglected phenomena (such as nonlinearities with 
conic sector bounds). This uncertainty is called “unstructured” because it is represented as 
a norm-bounded perturbation with no particular assumed structure. Plant uncertainty is 
called “structured” when there is real parameter uncertainty in the plant model, or when 
there are multiple unstructured uncertainties occurring at various points within the system 
simultaneously. Plant parameter uncertainty can arise from modeling errors (which usually 
result from assumptions and simplifications made during the modeling process and /or from the 
unavailability of dynamic data on which the model is based) or from parameter variations that 
occur during system operation. 

Robust control design and analysis methods for systems with unstructured uncertainty 
are accomplished via singular value techniques (refs. 1 through 9 and 21). For systems with 
structured plant uncertainty, however, a technique which takes advantage of that structure, such 
as the structured singular value (SSV) (refs. 16 through 27) or multivariable stability margin 
(MSM) (refs. 28 through 33), should be used. In order to compute the SSV or MSM, it is required 
that the system be represented in terms of an M-A model. The M represents a transfer function 
matrix M(s) of the nominal closed-loop system, and A represents an uncertainty matrix acting 



on M(s). The system M(s) is formed by closing the feedback system K(s) around the nominal 
open- loop plant interconnection structure P(s), as shown in figure 1. 



M(s) 

Figure 1. Block diagram of general M-A model. 


For multiple unstructured uncertainties, A is a block diagonal matrix, and for real parameter 
uncertainties, the A matrix is diagonal. As indicated in the literature (refs. 17, 18, 20, and 21), 
this model can always be formed for any linear interconnection of inputs, outputs, transfer 
functions, parameter variations, and perturbations. However, very little of the literature 
discusses methods for obtaining an M-A model. While formulation of an M-A model for 
unstructured uncertainties does not pose a major problem, forming an M-A model for real 
parameter variations can be very difficult. In reference 29, De Gaston and Safonov present an 
M-A model for a third-order transfer function with uncertainty in the location of its two real 
nonzero poles and in its gain factor. Although the given M-A model is easily obtained for 
this simple example, other examples do not yield such a straightforward result. A general state 
model of M(s) for additive real perturbations in the system A matrix (where A is assumed to 
be closed loop) is discussed in reference 34. Unfortunately, this model is not general enough 
for many examples, since system uncertainty is restricted to the A matrix and the uncertainty 
class is restricted to be linear. Morton and McAfoos (ref. 26) present a general method for 
obtaining an M-A model for linear (affine) real perturbations in the system matrices (A, B, 
C, and D) of the open-loop plant state model. In this model, the interconnection matrix P(s) 
is constructed first for separating the uncertainties from the nominal plant, and then M(s) is 
formed by closing the feedback loop. The M-A model thus formed can be used in performing 
robustness analysis of a previously determined control system. Alternatively, if the feedback 
loop is not closed, //-synthesis techniques (refs. 19, 20, and 21) can be applied to the P- A model 
for robust control system design of K(s). The result of Morton and McAfoos essentially reduces 
to that of reference 34 when the perturbations occur only in the A matrix (and the A matrix of 
ref. 34 is assumed to be open loop). An algorithm for easily computing M(s) based on the result 
of Morton and McAfoos is presented in reference 35. Although this method of constructing an 
M-A model is adequate for linear uncertainties, many realistic problems require a more general 
class of uncertainties, since, for real problems, uncertainties can arise in a nonlinear functional 
form (such as nth-order terms and cross terms). For these cases, it often becomes necessary to 
have repeated uncertain parameters in the A matrix, as is discussed later. Since the M-A model 
is a nonunique representation, models of various dimensions (due to the repeated parameters) 
can be constructed depending on how the model is obtained. It is therefore desirable to obtain 
one of minimal dimension so that the complexity of the SSV or MSM computations during 
robust control system design or analysis can be minimized. However, none of the literature (to 
the authors’ knowledge) addresses the issue of minimality. 
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This paper presents a methodology for constructing a minimal P-A model for single-input, 
single-output (SISO) systems with real parametric multilinear uncertainties, where the terms 
“minimal” and “multilinear” are defined as follows: 

Definition: A P-A model is minimal if the A matrix is of minimal dimension. It is 

shown later that A is minimal when its dimension is as close to the num- 
ber of independent uncertain parameters as possible, that is, when it con- 
tains a minimal number of repeated uncertain parameters (if any). 

Definition: A function is multilinear if the functional form is linear (affine) when 

any variable is allowed to vary while the others remain fixed; for exam- 
ple, /(a, b,c) = a + ab + be + abc is a multilinear function. 

The requirement that the P-A model be minimal provides a means of improving the efficiency 
of the SSV and/or MSM computations during robust control system design or analysis. The 
allowance of multilinear functions of the uncertain parameters provides a means of handling 
certain nonlinear terms in the transfer function coefficients (and, hence, the system matrix 
elements), namely, cross-product terms. In addition, nonlinear nth-order terms can be 
approximated within the multilinear framework, although this representation is conservative. 
The proposed procedure determines the minimal A matrix and the state-space form of P(s) 
given the system transfer function in terms of the uncertain parameters, where any or all 
the numerator and/ or denominator coefficients can be multilinear functions of the uncertain 
parameters except the leading denominator coefficient. It should be noted, however, that this 
procedure is presented without formal proof. Moreover, the state-space form used in modeling 
P(s) is an extension of the result of Morton and McAfoos for real parametric linear (affine) 
uncertainties (ref. 26). An extension of this result to multiple-input, multiple-output (MIMO) 
systems appears possible and is under study. The paper is divided into the following sections: a 
formal statement of the problem to be solved is presented in section 2, followed by a discussion 
of minimality considerations in section 3; the solution structure is presented in section 4, 
and computational details of the solution are presented in section 5; finally, several examples 
illustrating the solution are given next in section 6, followed by some concluding remarks in 
section 7. 

2. Problem Statement 

Given the transfer function of an uncertain system G(s, 6) as a function of the real uncertain 
parameters 6, find a minimal P-A model of the form depicted in figure 2 such that 

1. The diagonal uncertainty matrix A is of minimal dimension 

2. The model of the nominal open-loop plant P(s) is in state-space form 
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Figure 2. Block diagram of general P-A model. 
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The model must handle multilinear uncertainty functions in any or all the transfer function 
coefficients except the leading denominator coefficient, in order to construct a minimal P-A 
model, the dimension of the A matrix must be minimized. Hence, factors which have been 
found to affect the dimension of the P-A model are discussed next, followed by the approach 
used in forming a solution to this problem. 

3. Minimality Considerations 

In constructing a P-A model of an uncertain system, the A matrix can become unnecessarily 
large due to repeated uncertain parameters on its main diagonal, ft is therefore of interest to 
examine the factors which can cause this repetition, so that the number of repeated uncertain 
parameters can be minimized. A factor which can be shown to increase the size of the A 
matrix is the particular state-space realization used in representing the system. To illustrate 
this, consider the following simple examples: 


Example 3.1: Consider the system 

1 = 1 

^ T { 6 \ + $2) s + 0102 ( s d- ^l)( s + ^ 2 ) 


where 

01 = 01 o + h 02 = 02 o + 02 

-1 < <5i < 1 - l < h < 1 

and the system is represented in state-space form as 


A = 


0 

~ 0 \ 0'2 


1 

(-01 - 02 ) 



C = [1 0] 


This system has two uncertain real poles, 9 \ and 62 • The terms and 62 represent the 
uncertain parameters associated with the uncertain poles. In block diagram form, this system 
can be represented directly from the transfer equation 

s 2 y = -(0i + 02>y “ 0102Y + u 

= [-(01o + 02 o) “ (01 + 02)] S Y + (-01o02o - 02o01 ~ 01o02 - 0^2) y + u 


as shown in sketch A. 



Sketch A 
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Then the uncertain parameters can be separated out into the P-A model shown in sketch B. 



P(i) 

Sketch B 


The A matrix associated with this model is A = diag \8\, ^ 2 ]) which is four-dimensional. 

The nominal system P (5) can be represented as 
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y = [1 0] 
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Example 3.2: The system G(s ) in example 3.1 can also be represented by the block 
diagram in sketch C. 
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(s + & ] ) (s + 
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Sketch C 


5 








where 


#1 ‘ #1 o + ^1 

-1 < < 1 


#2 — @20 + ^2 
- 1 < 6 2 < 1 


and the state-space representation is given by 


A = 





C = [0 1] 


An equivalent block-diagram realization for this system is given directly by the cascaded form 
in sketch D. 



Sketch D 


An equivalent P-A model is given in sketch E. 

A 



P(s) 

Sketch E 


where the A matrix is given by A = diag which is two-dimensional, and a state-space 

realization of the nominal plant P(s) is given by 
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The P A models obtained in examples 3.1 and 3.2 are very different, even though they 
both represent the same system with two uncertain parameters. Obviously, the P-A model of 
example 3.2 is minimal (A is two-dimensional) and that of example 3.1 is not minimal (A is 
four-dimensional), since for minimality the dimension of A should be as close to the number 
of independent uncertain parameters as possible. These examples therefore demonstrate the 
effect the particular realization has in forming a P-A model. In particular, it appears that a 
cascade realization is a desirable form for obtaining a minimal P-A model. Thus, a general 
cascade-form realization is part of the approach taken in constructing a minimal P-A model 
so that uncertain real poles and zeros can be cascaded whenever possible. A problem arises, 
however, in that some transfer functions have a form which precludes cascading uncertain real 
poles or zeros such as 


G(s,6) 


frjS 2 + b^s + 63 
(s + #i)(s + 62 ) 


G(s,5) = 

s 1 + ais + a 2 


where 9 X and 0 2 are assumed to be uncertain (and hence a function of 6). Cascading the poles 
and zeros for either case would result in improper transfer function blocks being realized. For 
these cases, it is unavoidable for the minimal A matrix to have repeated uncertain parameters 
on the mam diagonal. However, for each inseparable pole or zero pair it is only necessary to 
repeat one uncertain parameter. This issue is addressed in the proposed solution, and a minimal 
P-A model for the first transfer function above is given as an example. 

Another factor which affects the dimension of the P-A model is the form of the coefficient m 
the system transfer function. If any of the coefficients are nonlinear functions of the uncertain 
parameters instead of multilinear functions (e.g., there are nth-order uncertain terms in any 
of the coefficients), then extra dependent uncertain parameters must be defined in order to 
represent these terms in a multilinear form. For example, Sj would be represented as 6i6 2 , 
where S 2 = <5f , and both ^ and S 2 would appear in the A matrix. Thus, for this case, it is 
again necessary that the minimal A matrix contain repeated uncertain parameters on its main 
diagonal. An example illustrating this situation is presented later. 

These issues are addressed in the solution presented herein for constructing a minimal P-A 
model. The approach taken in forming this solution is described in the next section. 


4. Solution Structure 

Based on the problem definition and the minimality considerations outlined previously, 
several issues are addressed in forming a solution to the problem of constructing a minimal 
P-A model given the transfer function of an uncertain system. First, a general cascade-form 
realization is found which can be used to obtain a minimal P-A model. Second, the minimal 
A matrix is determined for any uncertain system such that extra dependent parameters are 
assigned to account for inseparable pairs of uncertain real poles or zeros as well as nonmultilinear 
(e.g., squared) terms. (As stated previously, it should be noted that the representation of nth- 
order terms in a multilinear form is a conservative approximation.) Third, a method of obtaining 
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a state-space realization of P(s) for any uncertain system is found. Therefore, the proposed 
approach for constructing a minimal P-A model is given as follows: 

1. Obtain a cascade-form realization of the system so that the state-space uncertain model can 
be written as 

x = Ax + Bu 1 ^ 

y = Cx + Du J 


A - Aq + A A B = Bq + B A ) 

C = C 0 + C A D = Do + D A f 

The terms with the subscript o (A a , B 0 , C a , and D a ) represent the nominal matrix 
components, and the terms with the subscript A (A^, B A , C A , and D A ) represent the 
uncertain matrix components. 

2. Determine the minimal uncertainty matrix A as described in the problem definition and 
depicted in figure 2, where A is defined as follows: 

A = diag [<$! , (5 2 , 63, • ■ • , <5m] = diag [<5j, S D } = diag {6} (3) 


where _ 

A g R mxm <5/ G R mi 8 D g R mD 8 e R 


The vector of uncertain parameters is 

6 = [61, 


the partition of 6 containing independent parameters is 

8j = [<5i, <52, ^3) ■ • • , 8 mj \ 


and the partition of 8 containing dependent parameters is 

8 R = [^TU/ + l ’ ^771/ +2 > ^?7l/-t-3 ’ ■ ' ' i 8m £) 


where 

m 

771 / 

m D 

Also 


minimal number of uncertain parameters 

number of independent parameters given in G(s , 6) 

minimal number of dependent (or repeated) parameters 

P = Aq 


(4) 


where 

p uncertain parameters input to P(s), p G R m 
q uncertain variables output from P(s), q G R m 


Thus, since m/ is given and fixed, the minimal A matrix results when m D (i.e., the number 
of dependent (or repeated) parameters in the A matrix) is minimal (or zero, if possible). 


rm 
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3. Determine the state-space model of the nominal open-loop plant P(s) having the following 


q 

y 


x — Aqx + [B xp 

Bo] | 

[uj 


Cqx 

X + 

P>qp 

B)qu 


C o 

. D 2/P 

D 0 


where B X p, Cqx> D qp, D qu , and D yp are constant matrices, 
in the equivalent shorthand notation defined as follows: 





L U J ) 

Thus, P(s) can also be written 


p ll (s) 

P l 2 ( 5 ) " 


Ao 

B xp 

Bo 

P 2l( s ) 

P 22<»_ 


C qx 

. C o 

d 9P 

[ Dyp 

D qu 
Do . 


P ll( fi ) = = C qx (sI-A 0 )- l B xp + T)q p ' 

vfs) 

P 12( s ) = = Cq X (sI - A 0 ) l B 0 + T> qu 

q f ^ 

^21 ( 5 ) — u(s) = ~ Ao) *B X p + D yp 

P 22 to = = C 0 (sl - A 0 )- l B 0 + D 0 j 


( 6 ) 


(7) 


and with the notation of equation (6), the individual transfer function matrices P„ ( s ) can 
be expressed as 13 


p n« 


Ao 

Ba :p 

Cgx 

P>qp . 


p 2 i(«) 


Ao 

B xp 

Co\ 



Pn(s) = 


P 22( s ) = 


Ao 

Bo 1 


Cgrx 

g. 

0 



Ao 

B 0 1 

Co 

D 0 \ 


(8) 


This notation is used in this paper (and in the literature) to conveniently represent a 
transfer function matrix in terms of the state-space matrices of the system realization. 
The last term in equation (6) should not be viewed as a partitioned constant matrix 
but as the partitioned transfer function matrices defined by equations (7) and (8). This 
distinction is made in the notation through the use of solid (as opposed to dashed) 
partitioning lines. y 


It should be noted that this model is an extension of the result of Morton and McAfoos 
(ref. 26), where the T> qp matrix was required to be zero. In this paper, however, D qp is 
allowed to be nonzero in order to model the multilinear (cross-product) uncertain terms. 


The results for constructing a minimal P-A model via this approach are presented in the 
next section. 


5. Computational Details of Solution 

The proposed solution is presented in four parts. The results for obtaining a cascade-form 
realization of the uncertain system are summarized first, followed by the results for obtaining 
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a minimal A matrix and a state-space realization of P(s). Then, a summary of the overall 
procedure is presented. 


5.1. Cascade- Form Realization 

Given the transfer function of an uncertain system in terms of its uncertain parameters 
G(s 6) it is desired to realize the system in a cascade form of first- and second-order subsystems. 
Thus if the transfer function is given in unfactored form, the numerator and denominator 
polynomials must be factored into first- and second-order subsystems, where the second-order 
terms are only used to represent complex conjugate and inseparable pole/zero pairs. The given 
transfer function can then be represented as follows: 

G(s, 6) = K y {6 ) G c (s, 6) G r (s , 6) K u {6) (9) 


where K u and K y represent input and output gain terms, respectively, and G r and Gq represent 
the first-order and second-order transfer function components, respectively. Then 


G r (s , 6) = G Rk {s, 6) G RkX (s, 6)... G R2 (s, 6) G Rl {s, 6) 
Gq(s, 6) = G(j ( (s, S) Gc e _ l (s,S)...Gc 2 (s,6 ) G Cl {s,6) 

hi- I s + hi 


G Ri {8,6) = 


S + Otj 


„ , cs_ b 3i-2 s2 + ft 3i-l s + b 3z 

s 2 + a 2i -\s + a 2i 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 


and 


k number of first-order blocks 

£ number of second-order blocks 

Any or all these transfer function coefficients may be uncertain. The uncertainty may arise 
from either the coefficient itself being uncertain or from the coefficient being a multilinear 
function of one or more uncertain variables. Therefore, for either case, any of the coefficients 
may be a function of 6. Furthermore, the uncertain variables may have either an additive, 
e = e 0 + <5e, or multiplicative, e = e 0 (l + ^e)> form. 

The following cascade- form state- space realization of this system is proposed. 


Ar 

0 

B rK u 

B c Cr 

Ac 

BcD/jAu 

Ky B C C R 

KyC C 

KyD(jT> R K u 


(14) 


where 


A R 


A R i 

B R 2 C Ri 

B f?3 D ^2 C l?l 

B i?4 D /?3 D #2 C fll 


0 


r A § 

b R 3 c R2 

b Ra B R 3 C R2 


B R k _ l B R k _ 2 • -- B R2 C R 1 B Rk-l^ R k-2 • • ■ 
B R k B R k -i ■ ■ • b R2 C Ri B Rk B h - 1 • • • B Rs C R2 


0 0 
0 0 
0 0 
0 0 


A R k - i 

B R k c R k _ l 



10 



B R = 


a R 1 

B R 3 B R 2 ' D Ri 

b R4 b R3 t> R2 d Ri 


BRk-^Rk _ 2 -D fl 2 D i?1 


(16) 


C i?- [( D fl fc D i?fcl ...D /?2 C i?1 ) (D^D i?fe l ...D i?3 C i?2 )...(D i?fc C ftA; _ 1 ) (c fl J] (17) 


D R= [ D i? fe D^_ 1 ...D^ 2 D /?1 


(18) 


The Aq, H(j, Cq, and Dq matrices have the exact same form as equations (15) through (18), 
except that the subscripts R and k are replaced by C and i, respectively. The submatrices are 
defined from equations (12) and (13) as follows: 


A R j = ~ a i 
C R t = hi ~ a ihi-l 
where the terms a and j3 are defined in equation (12), and 


D fl. = hi-\ J 


(19) 


A Cj ~ 


0 

~ a 2j 


1 


a 2j — 1 J 

C Cj = [(*3j - a 2jhj-2) (hj-i - ^- 1 ^- 2 )] 


B C- 


D Cj = hj-2 J 


(20) 


where the terms a and b are defined in equation (13). The state-space realizations 
B R., C Rv D^| and |a^c, B^c, C j lead to the ith first-order and j'tli 
second-order transfer matrices G R (s, 6) and G^.(s,^), respectively, where i = 1,2, . . . ,k, and 
J = 1,2,..., 

The resulting cascade-form realization of the uncertain system is therefore obtained from 
equation (14) as 


o 

cs 

< 

L. 

B — 

B r K u 

_ b c c r a c. 

U 

b C B R K u_ 


C = [ KyB c C R KyC c ] D = [ KyD c U R K u } 


> 


( 21 ) 


This model is a general cascade-form realization for any uncertain open-loop SISO transfer 
function. The model does not, however, handle nonmonic denominator polynomials with 
uncertain leading coefficients. This would result in fractional (i.e., rational) matrix elements 
in the realization with uncertain parameters in the denominator of these elements. For real 
uncertain poles or zeros, two factors determine whether the first-order (real) or second-order 
(complex) block form should be used. The first is the nature of the uncertainty associated with 
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these terms, and the second is the form of the transfer function. If the real pole or zero locations 
are the uncertain parameters and the transfer function form allows these poles or zeros to be 
separated out, then the real block form should be used. If the transfer function form does not 
allow this separation, then the complex block form must be used. Furthermore, if there is a 
pair of uncertain poles or zeros that cannot be cascaded, then the resulting minimum A matrix 
will have a repeated parameter on the main diagonal for each inseparable pole or zero pair. 
Alternatively, if the coefficients of the second-order polynomial associated with the real poles 
are the uncertain parameters, then the complex block form should be used. These cases are 
illustrated in section 6 of this paper. The formulation of the minimal A matrix is presented 
next. 


5.2. Minimal A Matrix 

In formulating the minimal P-A model, the minimal A matrix must be determined first. 
The minimal A matrix is defined as in equation (3) with 

m = m/ + m£) (22) 

where mj is the number of independent uncertain parameters and mjj is the number of 
dependent uncertain parameters that must be added. The independent uncertain parameters are 
those defined in G(s, 6). However, as discussed previously, the dependent uncertain parameters 
are those independent parameters that must be repeated due to nonmultilinear terms in the 
transfer function coefficients and/or pairs of uncertain real poles or zeros that cannot be 
cascaded. Thus, for A to be minimal, mp (the dimension of 6p) should be minimized. It 
can be shown that if the system transfer function is formed from a given minimal P-A model 
of an uncertain system, the coefficients of the numerator and denominator polynomials will be 
multilinear functions of the uncertain parameters. Unfortunately, the converse is not generally 
true because of the dependence of the P-A model on the realization used for the plant. However, 
if the general cascade-form realization posed in this paper is used, the multilinear form of the 
transfer function coefficients can be used to establish that m = mj (i.e., mj) = 0) unless there 
are real uncertain pairs of poles or zeros that cannot be cascaded. Furthermore, it can be shown 
that if the coefficients of all the factors of the numerator and denominator polynomials are 
multilinear functions, then the coefficients of the expanded polynomials will also be multilinear. 
However, if there are nonmultilinear uncertain terms in the transfer function, then dependent 
parameters must be defined (and included in A) to represent the nonmultilinear term in a 
multilinear form. Moreover, if the nonmultilinear term is of the form S n , then n — 1 dependent 
parameters must be defined. If there are pairs of uncertain real poles or zeros that cannot 
be cascaded, then one additional dependent parameter must be added for each pair, and the 
dependent parameter can be either of the uncertain real parameters in the pair. Therefore, the 
number m, as determined by these rules, is the minimal dimension of the A matrix for the 
uncertainty class considered in this paper. Once this minimal dimension is determined, the A 
matrix can be defined as a diagonal matrix, as in equation (3), with the specified uncertain 
parameters on the main diagonal. Examples which illustrate these cases are presented later in 
section 6. 

5.3. State-Space Realization of P(s) 

Once the cascade-form realization has been determined, the system can be modeled as in 
equations (1) and (2), where the elements of A A , B A , C A , and D A are known functions of 
the uncertain parameters. Since any nonmultilinear terms have been redefined in a multilinear 
form when the minimal A matrix is determined, these matrices are multilinear functions of 
the parameters. In order to obtain a state-space model for P(s) as defined in equations (5), 
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expressions for these uncertainty matrices must be determined in terms of the matrices B xp , 

Cqx, D qp , D qu , and D yp from the model. With equations (4) and (5), these expressions can 
be determined as follows: 

A a = B xp A (I- D qp A)~ 1 C qx = B xp {I - A Dgp)' 1 A C qx 

B A = B . xp A (I- Dqp A)~ l D qu = B xp (I - A Dgp)' 1 A Dg U 

= A (I — Dq p A) l C qx = D yp (I — A Dg P ) _1 A C qx 

d A = D yp A (I -Dqp A) l D qu = D yp {I - A D 9p ) -1 A Dg M 

The inverse term makes computation of D qp very difficult. Furthermore, the matrix inversion 
can cause A A , B A , C A , and D A to have fractional (i.e., rational) elements with uncertain 
parameters in the denominator, which is not allowed in the uncertainty class being considered. 
Thus, it is desirable to represent this factor in terms of its Neumann expansion (ref. 36): 

(I — A Dqp ) = I + A Dqp + (A Dg p ) 2 + (A Dg p ) 3 + . . . (24) 

where the latter form in equations (23) has been assumed. Then equations (23) can be rewritten 
as 

a a = B xp A Cg X + B xp |a Dg P + (A Dg P ) 2 + (A Dg P ) 3 + . . . j A Cqx 

Ba = B xp A Dg M + B xp | A Dg P + (A Dg p ) 2 + (A D 9p ) 3 + . . .J A D qu 

C A = A Cqx + Dyp |a Dqp + (A D gp ) 2 + (A D ?p ) 3 + . . .j A C qx 

d A = D 2 yp A D gu + D pp | A D qp + (A D gp ) 2 + (A D gp ) 3 + . . .j A D qu 

The second group of terms add in the cross terms of the multilinear uncertainty functions. Each 
term in the series adds a higher order cross-product term. Since A A , B A , C A , and D A are 
multilinear functions with a finite number of terms, the Dg p matrix can be defined to have a 
special structure such that the infinite series of equation (24) can be replaced by a finite series. 
Hence, convergence of the infinite series of equation (24) need not be considered. This special 
structure for Dg p is given as follows: 

(A D 9p ) r+1 = 0 (26 ) 

Thus, 

(/ - A Dgp)- 1 = / + A Dgp + (A Dgp) 2 + . . . (A Dqp) r (27) 

where r is the order of the highest cross term occurring in A A , B A , C A , and D A ; that is, 

r = max (O a , O b , O c , 0 D ) (28) 

and O a , O b , O c , and O b represent the order of the highest order cross-product term in A A , 
B A, C A , and D A , respectively. That is, for a general uncertain m x n matrix M, 

0 M = max [order (m^); for all* = 1, 2, . . . , m and j = 1, 2, . . . , n] (29a) 
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where the order of each is the order of its highest order cross-product term, and cross- 
product term order is defined as 

order (<5i 62 63 • • • $i) — * ~ 1 (* — 1)2,..., m) (29b) 

Thus, the maximum value of r is rmax = m — 1, where m is the dimension of the A matrix. 
The required structure for to satisfy equations (26) and (27) is given as follows: 


1. da = 0 (i = 1,2, ...,m) 

2. If dij ^ 0, then for i = 1, 2, . . . , m and j = l,2,...,m: 


a. d 


■ji 


b- d i©lJ®l = 0 or rfj®2J©2 


= 0 or di 


i©(m— l),j©(m— 1) 


(30) 


where 0 represents modulo m addition (ref. 37) over the set {1,2,..., m}; that is, 


a 0 b — 

The desired equations can therefore be written as 


a + b ( a + b<m ) 

a + b- m (a + b > m) 


} 


A Cqx + B xp j 

A D qp + (A D qp)^ + . 

. (A DqpY | A Cqx 

A T)q U + B xp i 

A D qp + (A D qp) + . 

* (A Dg»p) r | A D qu 

A Cqx + D yp j 

A Dqp + (A n qp ) 2 + . 

. . (A Dgp) r j A Cqx 

A D q U + T)yp * 

A D qp + (A T) qp ) + . 

. . (A Dgp) r | A T) qu ^ 


(31) 


where r is defined by equation (28). Since the A^, B^, C^, D^, and A matrices are known for 
the given system, equations (31) can be used to determine B X p, Cqx, D q U , ^ yp- and Dqrp. Once 
these matrices are obtained, the state-space model of P(s) is determined. Hence, a minimal 
M-A model can also be formed. 

A procedure which summarizes the necessary steps in obtaining a minimal P-A model using 
these results is presented next. 

5.4. Summary of Procedure 

The following is a summary of the procedure implied by the preceding proposed approach 
for forming a minimal P-A model of a given uncertain system: 

1. Obtain the system transfer function in factored form. The coefficients of each factor should 
be a multilinear function of the uncertain parameters. If necessary, define new dependent 
parameters to represent any nonmultilinear terms in a multilinear form. 

2. Define the number of parameters in the A matrix, m, using equation (22). In so doing, 
determine if any new parameters are required to model inseparable uncertain real pole or 
zero pairs. If there are inseparable real pairs, either uncertain parameter in the pair must 
be repeated. 
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3. Define the minimal A matrix as in equation (3), using the independent parameters defined 
in the given transfer function as well as those defined in steps 1 and 2. 

4. Obtain a cascade-form realization for the system as a function of the uncertain parameters. 

5. Express the system matrices as in equations (2). 

6. Determine the maximum order of cross-product terms r in A A , B A , C A , and D A 
as defined by equations (28) and (29). Then A A , B A , C A , and D A have the form 
represented in equations (31), where D qp has the special structure of equation (30) required 
by equation (26). 

7. Express A A , B A , C A , and D A as 

A A = A A0 + A A1 + A A2 + • • • + A Ar 1 
B A = B A0 + B A1 + B A2 + • • • + B Ar I 
Ca = C A0 + C A j + C A2 + ■ • • + C Ar | 

Da = D A0 + D A1 + d A 2 + • • . + D Ar J 

where the subscript i represents the cross terms of ith order in each uncertainty matrix. 

8. The B X p, C qx , D pp , and T) q u matrices are found with the expansion described in reference 26 
for the uncertainty matrices having zero-order cross-product terms; that is, define 


M = 

a ao 

1 

o 

< 

CQ 


_C A0 

d ao_ 


— M] + M 2 <5 2 + • • ■ + 


(33) 


where the Mj matrices are appropriately partitioned. For the case of repeated parameters 
(due to inseparable real poles or zeros, or due to nonmultilinear functions), the 
matrix associated with the repeated parameter must be nonzero. These matrices can 
be decomposed into the product of appropriately partitioned column and row matrices as 
follows: 



i 

S 

to 

----- 1 


I 1 

§ 

II 

71 

iS 


1 

2 

~2 

tS>. 

i— 


(34) 


where M forms the ith column of B xp , M^) forms the ith column of T>yp, forms 
the ith row of C qx , and forms the ith row of D qu . Thus, 


'xp — [ 


M B2 ... 

■ M bJ 

y yp = I 

M Di, 

m d 12 ■ 

• M D lm ] 

KJ 

II 

2 

m£ 2 .. 

1 T 

■ M 0»] 

II 

2 

bb 

to 

m S 22 

••• M D 2m ] 


(35) 


9. Use the higher order cross terms of A A , B A , C A , and D A , as in equations (32), to 
determine the elements of the T) qp matrix. An augmented matrix equation can be formed 
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with equations (31). Begin with the first-order terms and specify as many elements as 
possible. Continue with the second-order terms, and proceed until all elements of T> qp are 
specified. Check T> qp to ensure that the required special structure of equation (30) and, 
hence, equation (26) is satisfied. 

10. Form the minimal P-A model as given in equations (3), (5), and (6) and depicted in figure 2. 
If the M-A model is desired, the feedback control system K(s) can be closed as discussed 
previously. 

It should be noted that the matrices M#., Mp., and M^., obtained in decomposing 

the matrices in equation (34), are not necessarily unique. A method of formalizing this 
decomposition for computer implementation is not addressed in this paper. However, an 
algorithm is presented in reference 35 which accomplishes this decomposition as an extension 
to reference 36. Some examples are given in the next section to illustrate these results. 

6. Examples 

The following examples illustrate the proposed procedure presented in section 5.4 for the 
various cases discussed in the preceding sections. 


Example 6.1: This example illustrates the construction of a minimal P-A model 

for an uncertain system whose transfer function contains an uncertain coefficient which is 
a multilinear function of the uncertain parameters as well as a coefficient which is a 
nonmult ilinear function of an uncertain parameter: 


0 ( 3 , 6 ) 


1 

s 2 + 2(ljs + uj ^ 


where 


C = Co + w — u o + 


This is a second-order system with uncertain complex poles. The uncertainty appears in the 
damping and frequency characteristics of the complex poles. The procedure in section 5.4 is 
used to determine a minimal P-A model as follows: 


n 

1. As given, G(s, 6) is in factored form. However, the constant coefficient u) is not a multilinear 
function of the uncertain parameters. Substituting for u in the above transfer function 
therefore yields the problematic term In order to represent this equation in multilinear 
form, the following dependent variable is defined: 

<$3 = 


so that 



2. There are two independent parameters, 6^ and <5^, and one dependent parameter, £3, due 
to the nonmultilinear term. Thus, from equation (22), 


m = mj + rnu = 2 + 1 = 3 
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3. The minimal A matrix can be defined as follows: 

A = diag[<^, <%] 

although ordering of the uncertain parameters in A is arbitrary. 

4. The cascade realization using equations (13) and (20) is determined to be 



Thus, it can be seen that for this example, the uncertain parameters occur only in the A 
matrix. 


5. Separate the nominal and uncertain parts of the system matrices as in equation (2): 


A = Ao + A a = u 2 1 + / u o 

L u o [ (-2a ) 0 8 u) -f 8^8%) -2(uj 0 8^ -f + 8^8^) 

B = B 0 + B a= [jl + [»' 


C — C 0 + C 'a =1 oj + 0 0 


6. Find the maximum cross-term order, r: 

°A = order (<5 W <5 3 ) = order (^<5^) = 1 
°B = 0 C = 0 D = 0 

r = max (1, 0, 0, 0) = 1 =» r = 1 

7. Expand as in equations (32): 

A AO Al -2uj 0 6u -2(u) 0 8^ + (, 0 6uj) _ + <5 W <5 3 -2 \6^8 W 

8. Find B xp , C qx , D yp , and D qu by using equations (33)-(35): 

A --°- J - --- = M 1( 5 c + M 2 <5u, + M 3 <5 3 

C A0 j D A0. 
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where 



0 

0 

1 

t 

1 

0 

Mi = 

0 - 

2 u>o 

I 

0 


0 

0 

T 

1 

1 

0 


0 

0 

1 

1 

1 

0 

m 2 = 

-2lv 0 

0 

t 

I 

0 


0 

_ 

0 

T 

1 

1 

0 


0 

0 

1 

1 

1 

0 

M 3 = 

0 - 

-2Co 

1 

I 

_J. 

0 


0 

0 

1 

1 

1 

0 


It is noted that M 2 and M 3 are the coefficient matrices associated with the | repeated parameter 
8 As indicated in the procedure, neither of these matrices can be a null matr . _> 

nonzero 'elements of the coefficient matrix associated with ^ have been split columnwise to 

form M2 and M3. 


Note: 


This is equivalent to reassigning the repeated parameters in A^q, B A0> 
Cah and D\n, to be columnwise independent prior to determining the 
Mi matrices. For this example, the reassignment can be accomplished as 


Follows: 


a A0 = 


0 

— 2 ujq8(jj 


0 


-2 {i0 0 S(^ + Co^3) 


order to be consistent with the M2 and M3 matrices determined 


above. 


Then 



0 u>o 


0 


m 2 



u>o 0 


0 


M 3 


il 


0 


0 Co 


0 
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Thus 


B xp — 


0 0 0 
-2 -2 -2 


— 


0 0 0 


Cgx — 


o 0; o 

uj 0 0 

0 Co 


— 


9. Find T) q p by using equations (31) and (32). Since r = 1 for this example, equations (31) can 
be written as 

^A1 = D qp] A C q X 

®A1 = B^fA D qp] A D q U 
CAl “ D yp[ ^ D q p ] A Cq X 
^A1 = D^[A Dgp] A Dg W 

These equations can be combined into an augmented matrix equation as follows: 


^AJ. j ®Al 

c ai | d A i 


B 


xp 


D 


yp 


A D 


qp 


CqX J Dg^ 


Substituting into this equation for A Al , B Al , C Al , D A1 , B xp , D yp , C qx , D qu , and A 
and solving for D qp yields 

0 0 O' 

0 


D QP ~ 


rf- 0 


0 


-1 

^uJq 


0 


which is consistent with the special structure required by equation (30). To see this, consider 
equation (30). As required by condition 1, d u = 0, d 22 = 0, and d 33 = 0. Now consider 
condition 2. First, d 2l ^ 0 requires that d u = 0 and that d 32 = 0 or d 13 = 0. Since d n = 0 
and d\ 3 = 0, this condition is satisfied. Similarly, d 32 ^ 0 requires that d 23 = 0 and that 

dl3 = 0 or d 2 \ = 0. Since d 23 = 0 and d\ 3 = 0, the condition is satisfied. Hence, the special 
structure of equation (30) is satisfied. 

10. The minimal P-A model as shown in figure 2 can now be constructed by using equations (3), 
(5), and (6). The realization of P (s) for the resulting model can be depicted as follows: 


X 1 


0 1 


x l" 


o 

o 

o 


Pi 


roi 

. x 2 . 


—Wq —2£ 0 uj 0 


[ x 2 J 

+ 

-2 -2 -2 


P2 

+ 

i 






k- 


.P3_ 


L J 


u 


Ql 

<12 

<13 


0 <jJq 

u 0 0 

Co 


0 


X 1 

x 2 


+ 


0 

1 


0 

0 


0 si 0 


pi 

P2 

P3 


y = [i o] 


x i 
. x 2 . 
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The A matrix is given by 




Example 6.2: This example illustrates the construction of a minimal P-A model for an 
uncertain system whose transfer function contains an inseparable uncertain real pole pair: 


G(a,S) 


bis 2 + b 2 s + ^3 

(s + 0l)(s + 0 2 ) 


where 


h = h 0 + fy>i 


b 2 = b 2o + S b2 b 3 = ho + S b3 


Q\ = 8i 0 + 6qi 9 2 = 9 2o + 6q 2 


1 Since the numerator is second order with uncertain coefficients, the uncertain real poles 
" in the denominator cannot be separated into the real cascade form. The denominator must 
therefore be expanded, and the second-order (complex) block must be used in the realization; 

that is, 

_ bis 2 + &2S + b 3 
U{8 ' 0) - s 2 + (01+02)8 + 0i0 2 

2. There are five independent uncertain parameters in this system which must all appear in the 
A matrix, namely, 

6l, b 2 , 63, 0i, and 0 2 

In addition, since 9\ and 0 2 are an inseparable pole pair, either 0\ or 9 2 must be repeated 
in the A matrix. (It can be shown that if % or S e2 is not repeated, the T> qp matrix wi 
not have the required structure of equation (30) and the higher order cross terms will not 
be modeled correctly.) Thus, from equation (22) 

m = 5 + 1 = 6 


3. The A matrix can be defined as 


A = diag 



fi 2 , % > b b2' S b3 


where 


s 2 = Hi 


4. The cascade realization can be constructed by using equations (13) and (20) and is given as 
follows: 


A = 

C = 


0 1 

-0102 ~(h+H) 


B = 
D = 


0 

1 

6l 


(63 - 0 id 2 bi ) ( b 2 - { 9 1 + #2)^1) 

For this example, the uncertain parameters occur in the A, C, and D matrices. 
5. These system matrices can be expanded as in equations (2) to yield. 


20 


where 



0 1 
—0\o^2o ~(^lo + ^2o). 

[ (^3 o ~ ^lo^2o^lo) *2 — {*1 o + ^2o)^lo 


a a = 

C A = 


I 0 O' 

L ~ (^2o% + hoSto + *61*62) ~(*6l + *62) . 

C A11 C A12 


B 0 


0 

1 


D 0 = 


ho 



C A11 = *b3~*2oho*61~01oho*62^*lo^2o*bl~ho*61*62~01o*62*bl~^2o*61*bl~*ei*62*bl 

C A12 ~ *b2 ~ ho*6\ ~ *lo*82 ~ (^1 o + ®2o)*bl ~ ^ 01^61 — *62*bl 
6. Find the maximum cross- term order r using equations (28) and (29): 

O a = order {8q\8q 2) = 1 Oq = order (0) = 0 

Oq — order (<5^i^ 2 ^6l) = 2 Op = order (<S bl ) = 0 

r = max (1, 0, 2, 0) = 2 => r = 2 


7. Expand A A , B A , C A , and D A as in equations (32): 


where 


A A a A q + A Ai + A A2 B a - B ao + B A1 + B A2 

C A = C A0 + C A1 + C A2 d A = D A0 + D A1 + D A2 


A A0 

A A1 

A A2 


0 0 

. ~{® 2 o * 6 \ + ^10^2) ~(*61 + * 62 ) 

0 O' 

-*6i*62 0. 

0 O' 

0 0 


Bao - B A1 — b A2 = 


0 

0 


C A0 - [{*63 - ^lo^l - ho*lo*62 ~ hoho*bl) {*b2 ~ ho*81 ~ b lo*62 ~ Vl o + *2 o) *6il] 
C A1 = [~ { b lo*61*62 + 02o*6l*bl + ^lo^2^6l) - (*61*bl + *62*bl )] 


C A2 - [~*6l*62*bl °] 


°A0 = foil D a1 = [0] D A2 = [0] 
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8. Solve for B xp , C qx , D yp , and D yu using equations (33)-(35): 

M = Mifyi + M 2 <5 2 + M 3 6 92 + M 4^ + M b S b 2 + M 6^ 3 

where 



0 

1 

0 ! 

0 


- 0 ' 



Mi = 

-ho 

1 

1 

_°_J 

0 

= 

-1 


ho ojo 


_ ~^2o^lo 

I 

0 ! 

0 


- —b\o - 




1 

o 

o 

1 

1 

l 

ol 

r 

0 1 

r 

. i 


M 2 = 


0 b\ Q 


m 4 = 



0 


0 

i 

i 

i 

0 

m 3 = 

ho 

-1 

i 

i 

i 

0 


_ —8\oblo 

—ho 

i 

i 

i 

0 

0 

0 


1 

t 

1 

0 


0 

0 


1 

1 

1 

0 

== 

01o#2o 

-(^l 0 + ^2o) j 

i 




0 

1 

0 I 

0 


m 5 = 

0 

! 

_°J- 

0 





i 

i j 

0 

i 



0 

i 

0 ! 

0 


m 6 = 

0 

i 

0 

0 

= 



1 

n 

0 i 

0 



-1 
-ho J 

0 


0 1 ! 0 


-1 


ho ijo 


o 

1 


r 

1 

1 O 

l 


Oloho 

L -1 . 



1 

1 

0 o 1 

1 

1 


0 1 jo 

1 J 


- 


r o i 


_ 

1 

1 o 
1 


1 ojo 

i j 


- 


It is noted that Mi and M 2 are associated with the repeated parameter in A. As 
discussed in example 6.1, neither Mj nor M 2 can be a null matrix. Thus, as in example 6.1, 
Mi and M 2 were formed by splitting up the columns of the coefficient matrix associated 
with 6 ei . (It can be shown that if Mi and M 2 are not both required to be nonzero, B qp 
will not have the required structure, and hence the higher order cross-product terms will not 
be modeled correctly.) 

Note- As mentioned in the previous example, this is equivalent to reassigning the repeated 
parameters in A A0 , B A0 , C A0 , and D A0 to be columnwise independent prior to 
determining the M; matrices. For this example, the reassignment can be accomplished 
as follows: 
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a ao = 


I 0 0 

L ihohi + ^10^2) - ih + %) 


C A0 [t%! HMn e ^dn-ho«2ohi) {h2-hoh~hoie2-{ho + «2o)hi) 

in order to be consistent with the M, and M 2 matrices determined above. 

Thus, based on the above, the following results are obtained: 


n 

0 0 

0 0 0 

O' 

d 

II 

- 


*-*xp 

-i -1 

-100 

0 

J 

- ho - ho - ho - 1 1 1 


^2o 

0 


l 


j 


0 

1 



■ 0 ■ 


Cqx = 

e lo 

h o e 2o 

1 

{ 0 \ o + 0 2o ) 


^qu = 

0 

0 

-1 



0 

1 



0 



1 

0 



. 0 . 
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The higher order cross terms 
D qp as follows: 


are used as in equations ( 31 ) and ( 32 ) for r = 2 to determine 


A A1 + A A2 = Bxp [A D qp + (A D 9p ) 2 ] A C qx 
B A1 + B A2 = B xp [A D qp + (A D 9p ) 2 ] A D qu 
C A1 + C A2 = D pp [A D qp + (A D 9p ) 2 ] A C qx 
D A1 + D A2 = T> yp [A B qp + (A D 9p ) 2 ] A D qu 

Substituting for these known matrices and solving for U qp yields the following result: 

0 0 0 0 O' 

00000 
00000 
110 0 0 
00000 
00000. 

which satisfies the required structure of equation ( 30 ). 


D gp - 


’ 0 

0 

1 

$ 2 o 

1 

0 

0 
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, oN /cn ,i (a\ flip fnllowine P-A model can be formed as depicted 
10. Thus, from equations (3), (5), and (6), the following r u 

in figure 2 where the state-space form of P(s) is given y 

0 0 0 0 0 0" 

_1 -1-10 0 0 


x = 


0 1 
-01o#2o -(^lo + ^2o)J 


x + 


P + 


u 


02o 

0 


' o ' 


0 

1 


0 


^lo 

1 

X + 

0 

u + 

^lo^2o 

01o + #2o 


-l 


0 

1 


0 


1 

0 


_ 0 



0 

0 

1 

$2o 

1 

0 

0 


0 0 0 0 0' 
0 0 0 0 0 
ooooo 
110 0 0 
ooooo 
ooooo 


y = l^3o — ^lo^lo^lo 
+ [~Ho — ho 
and the A matrix is given by: 


62o-( 0 lo + 0 lo) b lo] x + M u 
-b l0 -11 1]P 


A = diag 


Sq \ , Sq \ i 6q 2 , 6 b \, ^62’ fy>3 


Example 6.3: This example illustrates the construction of a minimal P-A model for an 
uncertain system which has many independent uncertain parameters. 

(frs + /%) (03 s + &) ( hf^hs + h) 

G ( s ' 6 ) = (7+^j(I+ a 2 )(s 2 + ais + 02) 

where 

ai = aio + 6 al «2 = «2o + 6 a 2 «1 = a lo + S al a 2 = a 2o + S a2 

01=010 + 6(31 02 = 020 + 6(32 03 = 030 + 6(33 0A-0Ao + 6(3A 

b± = bi 0 + S bl b 2 = b 2o + 6 b 2 ^3 = b 3o + 6 b3 

1. The given transfer function is in correct form, all coefficients are multilinear, and all pole-aero 
pairs can be cascaded. 

follow: m = rnj + mo = 1 1 + 0 = 1 1 =>m — 11 

3. The minimal A matrix can be defined as 

A = diag [s a i, 6 a 2 , S ah S a2 , 6 01 ,_6 0 2 , 3 , ty 4 » 6 bl , S b2 , S b3 \ 


_ b i s 2 + b^s + 1?3 
^ _ s 2 + ais + a 2 


G R = G Rl G R2 


\0\s + 02] 

'03 s + 0A 

[ s + «i j 

L s + a 2 . 
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The realization of G(j is determined from equations (13) and (20) for ( = 1: 


A C = 

C C = 


0 1 

- 0-2 ~ a \ . 

63 - a 2 f>i 6 2 — a \b\ 


B c = 
D C = 


0 

1 


h 


and the realization of Gr is determined from equations (12) and (15)— (19) for k = 2: 


A R = 

CM 

1 

i ' 1 

1 1 

CM 

L- 1 L 


II 

QS 

' 1 ' 
A. 

Cr = 

I 03 (02 - a 10l) 04 - oi 2 03 

11 

q; 

Q 

0301 


Then the cascade realization is given by equations (14) and (21) for unity gains at the input 
and output: 


Thus 


A = 

a r 

B C Cr 

V 

O 0 

B = 

BR 1 
B C Dr\ 

C = 

D C Cr 


D = 

\d C Dr ] 


A = 

C = 


— oq 0 0 

02 ~ a\0\ -a 2 0 

0 0 0 

-03(02 ~ a 10l) 0A~ a 203 ~ a 2 

f h03(02 ~ a \0\) h(0A ~ <*203) 


0 ' 

0 

1 

~ a \ _ 

(63 - a 2 bi) 


(62-0161) 


1 1 



W3J 


D = 


0i03h 


5 9- For this example, uncertainty arises in the A, B, C, and D matrices. The A A , B A , C A , 
and D a matrices are fairly complicated for this example and are therefore not given. The 
order of the highest cross-product term is 3, so that r = 3. Following the procedure outlined 
in section 5, the results are determined in a straightforward manner to be 
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10 Equations (5) and (6) can now be used to obtain the state-space model of P(s) and the A 
matrix determined as step 3 of this example. This yields the desired P- A model as depicted 

in figure 2. 


These examples illustrate the proposed procedure for forming a minimal P-A model of an 
uncertain system. Although all the steps involved in obtaining these results have not been 
included (particularly in example 6.3), the stated results should provide a guide in performing 
the steps of the proposed procedure. It should be noted that, for ease of hand computation, 
examples 6.2 and 6.3 included only the simplistic (and less realistic) case in which the coefficients 
themselves are the uncertain parameters. However, it is emphasized that, as illustrated in 
example 6.1, the proposed procedure does handle the more realistic case in which the uncertain 
transfer function coefficients are multilinear functions of the uncertain parameters. 

7. Concluding Remarks 

A proposed procedure is presented for forming a P-A model of an uncertain system which 
appears to be of minimal dimension, given its transfer function in terms of the uncertain 
parameters. The uncertainty class considered in this paper allows the transfer function 
coefficients to be multilinear functions of the uncertain parameters, and the uncertainties may 
arise in any or all of the A, B, C, and D matrices of the system model. Although no proofs are 
presented regarding minimality, the resulting models appear to be minimal m dimension for all 
examples worked thus far. Moreover, even if some counterexample exists for which the resulting 
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P A model is not minimal, the outlined procedure does provide a means of handling the more 
realistic uncertainty class which includes multilinear functions of the uncertain parameters. This 
procedure involves realizing the system in a cascade form, determining the minimal A matrix 
of uncertain parameters, and obtaining a state-space model for the nominal open-loop system 
P(s). As stated previously, the minimal M-A model can then be easily obtained by closing 
the feedback loop. Three examples were given to illustrate the proposed procedure. The first 
example had uncertainty in the A matrix only and illustrated how this method handles uncertain 
system transfer functions having coefficients that are multilinear functions of the uncertain 
parameters as well as coefficients that are nonmultilinear functions of the uncertain parameters, 
for which a repeated uncertain parameter had to be used. This example is representative of 
typical problems that can arise in the mathematical modeling of realistic dynamic systems. 
The second example had uncertain parameters arising in the A, C, and D matrices only. This 
example illustrated the formulation of a minimal P- A model for a system with inseparable real 
uncertain poles, and also involved repeating an uncertain parameter in the A matrix. The third 
example had 11 independent uncertain parameters, which arose in the A, B, C, and D matrices 
of the system realization. This example was included to demonstrate this method for a large 
number of uncertain parameters which occur in all the system matrices. Also, for this example 
none of the uncertain parameters had to be repeated. Thus, the minimal A matrix contained 
only the independent uncertain parameters given in the problem. 

_ Further work on the proposed procedure could include systems having a nonmonic character- 
istic polynomial with an uncertain leading coefficient, as well as systems having an inner feedback 
loop which may or may not have uncertainties. The former case might require extending the un- 
certainty class to include rational expressions containing multilinear functions of the uncertain 
parameters in the numerator and denominator, and the latter case might require a modification 
in the formulation of the cascade realization. Although the procedure presented in this paper is 
for single-input, single-output systems, an extension to multiple-input, multiple-output systems 
appears possible and should primarily involve modifying the cascade- form realization. Other 
areas of future work include development of a means of proving the minimality of a given P-A 
model, development of a proof that the above procedure yields a minimal P-A model, and 
development of a method of reducing a nonminimal P-A model to a minimal form. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
May 10, 1991 
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